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We present a set of equations describing the evolution of the scalar-type cosmological perturbation 
in a gravity with general quadratic order curvature coupling terms. Equations are presented in a 
gauge ready form, thus are ready to implement various temporal gauge conditions depending on the 
problems. The Ricci-curvature square term leads to a fourth-order differential equation for describing 
the spacetime fluctuations in a spatially homogeneous and isotropic cosmological background. 

PACS numbers: 04.50. +h, 04.62. +v, 98.80.-k, 98.80.Hw 

1. Introduction: In a series of work jl],f| we have been studying the evolution of cosmological perturbations in a 
gravity with additional curvature-square correction terms 



S = 



AR 2 + BR ab R ab ) + L r , 



(1) 



where L m is the matter part Lagrangian. The gravitational field equation is presented in Eq. (2) of Q|. The one- 
loop order quantum correction generically leads to curvature-square terms, and due to Gauss-Bonnet theorem the 



In the literature, many works analyzed the roles of 
. However, we find rarely any work considering the 



two terms in Eq. (|lj) completely account such corrections, 
R 2 term in the context of cosmological structure formation, 
Ricci-curvature square term in the same context. When the curvature-square terms come from quantum correction 
of Einstein gravity, if R 2 term is important R ab R a b term should be important as well because two terms are of the 
same order. 

As a background world model we consider a spatially flat, homogeneous, and isotropic, Friedmann-Lemaitre- 
Robertson- Walker (FLRW), model. The general perturbations of the FLRW spacetime can be decomposed into 
three-types based on the tensorial properties of the perturbed variables: these are the scalar-type, (transverse) vector- 
type, and (transverse-tracefree) tensor-type perturbations. The tensor- and vector-type perturbations correspond to 
the gravitational wave and the rotational perturbation, respectively. The scalar-type perturbation includes the den- 
sity condensation and the corresponding longitudinal type velocity perturbation, gravitational potential (curvature) 
perturbation, etc. Due to the high symmetry in the background model (homogeneity and isotropy) the three-types 
of perturbations decouple from each other in the linear level and evolve independently. Thus, without losing general- 
ity we can consider the three perturbations independently. The gravitational wave and the rotational perturbations 
were analysed in |Q and Q , respectively. For the gravitational wave we derived a fourth-order differential equation, 
whereas, the rotational perturbation is simply described by the conservation of angular momentum. In this Brief 
Report, we will present a complete set of equations for the scalar-type perturbation in a gauge-ready form: these are 
in Eqs. (jj|]l4|). We will show that, in general, the Ricci-curvature square term leads to a fourth-order differential 
equation. Applications to specific cosmological scenarios are left for future work. 



2. Scalar-type perturbation: The cosmological spacetime metric with the most general scalar-type perturbation is 

ds 2 = -a 2 (1 + 2a) dri 2 - a 2 (3, a dr)dx a + a 2 [5 a0 (1 + 2(p) + 2% a0 ] dx a dx< 3 , (2) 

where a(r\) is a cosmic scale factor. a(x, 77), /3(x, 77), tp(x,rj) and 7(x, 77) are perturbed order metric variables. Except 
for two degrees of the gauge redundancy, these variables represent the scalar- type metric perturbations. It is convenient 
to introduce the following combinations \ = a (f3 + a i) and k = 3 (Ha — <p) — (A/a 2 )x', A is a Laplacian based on 
5 a p, an overdot denotes the time derivative based on t where dt = adrj, and H = a/ a. a is the perturbed parts of 
the comoving lapse function. tp, x^ an d k are the perturbed parts of the three-space curvature, the shear, and minus 
of the perturbed expansion scalar of the normal frame vector field, respectively; see Sec. 2.1.1 in We decompose 
R(x,t) — R(t) + SR(x,t); unless necessary, we ignore the overbar which indicates the quantity to the background 
order. For later use we present: 

R% = 3 (H + ff 2 ) - (k + 2Hk + 3Ha + ^a^j , R° a = (k+ 




1 



k + 6Hk + 3Ha + — {3(p + x + H X ) 



(3) 



where = rj and | indicates a covariant derivative based on 5 a p. 

3. Gauge strategy: The gauge transformation properties of the metric variables are presented in Sec. 2.2 of [pj. 
Under the gauge transformation x a = x a + £° with T = a£°, we have [see Eq. (17) in ||]: 



a-T, fi = <p- HT, x = X-T, 



n+ (3H + A/a 2 )T, SR = SR-RT. 



(4) 



The variables a, ip, X i K ) an d SR only depend on the temporal gauge transformation, thus are spatially gauge invariant; 
for 5R we used Eq. (^), see also Eq. (p|). We have several choices for the temporal gauge fixing: the synchronous 
gauge (a = 0), the uniform-curvature gauge (tp = 0), the zero-shear gauge (x = 0), the uniform-expansion gauge 
(k = 0), and the uniform- R gauge (SR = 0). Unless we consider perturbations in the energy- momentum tensor we do 
not have the comoving gauge (or the uniform- field gauge for the scalar field), the uniform-density gauge, etc., which 
are related to imposing conditions on ST£. Except for the synchronous gauge condition, each of other temporal gauge 
conditions completely fixes the temporal gauge mode. Thus, a variable in such a gauge condition corresponds to a 
unique gauge invariant combination involving the variable concerned and the variable used in the gauge condition. 
We proposed the following notation for the gauge invariant variables || : 



R R 

ip x = (p-H X =-H X ip, 6R X = 6R- Rx = -Rxsb., 5R V = SR — —ip = ——cpsn, etc. 



(•5) 



In this manner we can construct systematically all possible gauge invariant combinations. Since there exist several 
(in principle, infinitely many) different gauge conditions available for each variable, and correspondingly several gauge 
invariant combinations, our notation for the gauge invariant variable is convenient in practice. 

In the Ricci-curvature square (we call it Ricci-square) gravity we have several additional temporal gauge conditions. 
From Eqs. (^,^) we have: 



6 EX = 5R° n 



R U T, 



SR" = SR" 



-HT, 



5R<$ 



8R<$ 



1 



R^SpT. 



(G) 



Thus, a temporal gauge condition among 5R° = 0, SR 



0. 



and SR^ 



SR" — also completely fix the temporal 



gauge degree of freedom; but, these are combinations of our fundamental gauge conditions, because these conditions 
are the same as SR — A[Hn + (A/a 2 )ip] = 0, n + (A/a 2 )x = 0, and Hk + (A/a 2 )ip = 0, respectively, 

In the following section we will present the perturbed set of equations without imposing the temporal gauge condition. 
In this way, the right to impose the temporal gauge condition can be used as an advantage for handling problems 
depending on the mathematical simplification or the physical interpretation we can achieve Often, different 

gauge conditions suit for different problems and keeping the equations ready for imposing various gauge conditions 
(thus, in a gauge ready form) is practically convenient. 



4- Equations in the gauge ready form: Equations for the background are presented in Eq. (7) of fi]]: 



H 2 + 2 (3A + B) (2HH - H 2 + 6H 2 H ) = -/i, 



/' : 



-3H(n + p), H=-—[T-T + v+p), (7) 



where T = 1 + 2(A + B/S)R, \i = -T °, and p = \TJ; T ab is the energy-momentum tensor of the additional matter 
part Lagrangian. The third equation follows from the first two equations. 

The perturbed set of equations follows from the gravitational field equation in Eq. (2) of Q. The complete set of 
equations in the gauge ready form is the following [we introduce F = 1 + 2AR, thus 5F = 2ASR]: 
Definition of n: 



. A 
-3 (ip -Ha) -x- 



Energy constraint: 



ST$ = 2F 
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3IH + H' 
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7a 



SF 
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Bi2[HK+^r<p) +6H [Hk+ -ftp 
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2{H-6H*)-- S 
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a- 



A 



Hk+-z<p)+ 3HSR -[3H' i + -f)5R 



A' 



+6HHa + 6 (-HH + 2H 2 - 9F 2 ij) a - (&H + 5F 2 ) k - -H— + - 
Momentum constraint: 



A 

2* 



(9) 



T° - --V 

3 a 



_/ A 3F 



A 



(<5F - HSF^j 



-B{ - ! /• H- — A ) - :!// U+^X 



-3ifd- 



(2F + 9F 2 ) 



a + 3- 



A 



2#> + (ij - 2F 2 ) x 



(10) 



Tracefree propagation: 
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(11) 



Raychaudhuri equation: 



<5F^ - (5F ° = 2F < 
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2SR + 6H 2 SR 



3H 2 H 
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Trace equation: 



and we have 



ST = -SR -2{3A + B) 



5R 



A 

a 2 



SB + 3I.JSR - —SR -R( K + a)-(2R + 3HRJ a 



k + AHk + 3Ha + — (2tp + a) 
a 2 



(12) 



(13) 



(14) 



We decomposed the equation for 8TS into tracefree and trace parts. Equations (|^-[l4|) form a redundantly complete 
set describing the scalar-type perturbation in a gauge ready form; for example, the combination of Eqs. (|^Jl^) leads 
to Eq. (H). In the limit of R 2 gravity (thus 5 = 0) Eqs. (|[]l|) correspond to Eqs. (54-58,60) in @ which also 
correspond to Eqs. (22-26) in J5j; the names of the equations are based on our convention in J5j. When we analyse 
we have the right to impose one temporal gauge condition out of several choices displayed in Sec. 3. 

5. R 2 gravity: We set B = 0. Without additional matter terms (T a b = 0), the analysis can be made most 
conveniently in terms of ip in the uniform-F gauge where we take SR = (thus SF = 0) as the gauge condition; 



3 



equivalently, we let SR = and change the other perturbation variables into the gauge invariant forms, for example, 
ip — > (fisR, etc. Using Eq. (|8|-|l0|) we express k and a in terms of (p. From Eq. ( |l3|) we can derive 



(HF/F + 1/2) 2 



a 3 F 

;¥>sr 



{HF/F + 1/2) 2 



-~<PSR = 0. (15) 



This equation was derived in a more general context in ; it is valid for a general action with f(R) replacing R+AR 2 
with F = df/(dR). Perturbation analyses in the R 2 gravity were also pursued in j| in different gauge conditions. 
Using ipsR — —{H/R)5R V in Eq. (||) we can derive the equation for 5R V . In the large scale limit (thus, ignoring the 
Laplacian term) we have an integral form solution valid for the generally time varying background dynamics 



<p SR (x, t) = -:§ 6R v (x, t) = C(x) + |d(x) / ( " HF/F 3 t 1/2)2 dt, (16) 
R o Jq a r 



where C(x) and D(x) are integration constants indicating the growing and decaying solutions, respectively; the 
decaying solution is higher order in the large scale expansion, see Eq. (109) in 0. Thus, the non-transient solution 
of ifsR is conserved in the large scale limit, and the generalized nature of the gravity does not affect this result! 



6. Discussions: In general we expect the curvature square terms lead to fourth-order gravity theories. It happens 
that, as we just saw, the R 2 or f{R) gravity leads to a second-order system. This is probably due to the conformal 
symmetry of the f(R) gravity to the Einstein gravity with a minimally coupled scalar field, The Ricci-square 
term does not have such a conformal symmetry, and we expect to have a fourth-order differential equation for the 
scalar-type perturbation. In general, we do not know which gauge condition will be the most suitable for the problem 
a priori. As a trial attempt we can impose the same gauge condition used in R 2 gravity which is the uniform- R gauge; 
thus, we set 5R = and 5F = 0. One possible way to derive a fourth-order equation is the following. Using Eqs. 
(I§|, [l3||l4"l) we can express a, d, and k + (A/a 2 )x in terms of k and <p. From a time derivative of Eq. (|l4|), and from 
Eq. (Bfiwe can derive a coupled set of two second-order differential equations for k and Hk + (A/a 2 )ip. In such a 
special form, however, it happens that the equations alone cannot be reduced to Eq. ( |l5| ) in a pure R 2 gravity limit. 

Our main result is Eqs. (pfjl^) which describe the evolution of a perturbed flat FLRW background in a gauge 
ready form; the evolution of background is governed by Eq. (0). The set of equations for the Ricci-square term is 
a new result, and after imposing a suitable gauge condition it generally leads to a fourth-order differential equation 
without remaining gauge mode. We mentioned one possible set of a fourth-order differential equation in a previous 
paragraph as a demonstration. However, the mentioned set may not be the best form possible in the available pool of 
the variables and the gauge conditions in Sec. 3. Equations ([sHli]) are arranged in a general gauge-ready form, thus 
allow an easy implementation of one's favored choice for the variables and the temporal gauge. 

The fourth-order terms in Eq. (|l|) usually appear as the one-loop order quantum corrections. In this regard such 
terms can be regarded as a transient (because higher-order correction terms will appear soon in higher order quantum 
corrections) form of medium-energy effective action || . Thus, we regard the A and B terms in Eq. ([!]) as corrections 
to Einstein action, and in such a case, we should consider the third and fourth order derivatives as corrections to the 
second order equation based on Einstein gravity. In the case we should accept Eq. ([!]) as the fundamental theory 
we are not certain whether the possible four solutions of the expected fourth-order differential equation in a general 
situation (and also the fourth-order gravitational wave equation derived in (!]]) should be interpreted as representing 
physical degrees of freedom; for example, we should consider the Cauchy problem and the stability of the full theory. 
We will not address this important question in this Brief Report, ]l0| ] . 

We briefly describe how we can connect specific solutions of our formulation with the present day observable 
quantities. Since we have some conceptual difficulties mentioned in the previous paragraph for fourth-order gravity, 
let us consider the R 2 gravity as an example. Let us consider a scenario where the higher-order gravity dominates 
the early evolution stage of our observable patch of the universe, and at some point Einstein gravity takes over the 
dominance till the present era. If we further assume that the higher-order gravity era provides an accelerated expansion 
(inflation) stage, the observationally relevant scales may transit from subhorizon to superhorizon scales during the 
stage. If conditions (e.g., allowing an analytic form mode-function solution of the perturbed equation) are met, we 
can derive the generated quantum fluctuations (usually, in general scales) based on vacuum expectation value, and as 
the scale becomes superhorizon it can be interpreted as the classical fluctuations based on the spatial average. As long 
as the relevant scales remain in the superhorizon stage during the transit epoch of the gravities (from higher-order to 
Einstein) there exists a conserved quantity, ip in the comoving gauge is the conserved quantity in the superhorizon 
scale independently of the changing gravity (as long as the gravity belongs to a certain class), changing field potential 
(for the scalar field), and changing equation of state (for the hydrodynamic situation), see In our case, ip$n in 
Eq. ( |l~6| ) is the conserved quantity; the growing solution is conserved as C(x). Since we are considering the linear 
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perturbation all quantites are related to each other by linear combinations; for the specific relations to the quantities 
with observational consequences (e.g., temperature anisotropies in the cosmic microwave background radiation, and 
the Newtonian hydrodynamic variables) see Eq. (23) in B. In summary, the curvature fluctuation (p in a certain gauge 
is conserved and connects perturbations in both eras. In later Einstein era with hydrodynamic contents, the density 
fluctuations will fall-in to the potential well of, thus co-grow with, the curvature fluctuation. Concrete analyses of 
the quantum generation and classical evolution processes with the generated observational specta in various scenarios 



based on a class of generalized gravity theories were carried out in 1 1 1 . Applications of our formulation to specific 
cosmological models are left for the future work. 
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